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ABSTRACT 



in 

O The merger of two white dwarfs (a.k.a. double degenerate merger) has of- 

^ ten been cited as a potential progenitor of type la supernovae. Here we combine 

I ^1 population synthesis, merger and explosion models with radiation-hydrodynamics 

light-curve models to study the implications of such a progenitor scenario on the 
>■ observed type la supernova population. Our standard model, assuming double 

degenerate mergers do produce thermonuclear explosions, produces supernova 



^O light-curves that are broader than the observed type la sample. In addition, we 

O 

o 
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discuss how the shock breakout and spectral features of these double degenerate 
progenitors will differ from the canonical bare Chandrasekhar-massed explosion 
models. We conclude with a discussion of how one might reconcile these differ- 
ences with current observations. 

Subject headings: Nucleosynthesis, Stars: Supernovae: General 



Introduction 



Type la supernovae are among the strongest explosions observed by astronomers. These 
thermonuclear bombs dominate the production of many iron-peak elements in the universe. 
Their near-uniform and easily calibrated light-curves have allowed astronomers to use them 
as probes into the universe and the nature of cosmology. However, despite their extreme 
importance, we still do not understand the details of their explosion mechanism nor do 
we understand the exact nature of their progenitor. Although it has been argued that 
white dwarf systems with masses below the Chandrasekhar limit may produce thermonuclear 



explosions and, in WD- WD collisions, this is generally true (Raskin et al. 2009 Rosswog et 



aL]|2009 ), it is generally believed that most type la supernovae are produced as an accreting 
WD is pushed above the Chandrasekhar mass limit (see Livio 2001 for a review). In this 
paper, we study the progenitors of these Chandrasekhar-massed explosions. 

The exact nature of the binary progenitor that leads to the WD accretion is also a 
matter of debate. The wide range of type la progenitors have been divided into two clases 
(see Livio et al. 2000 for a review): double degenerate mergers where the white dwarf 
accretes material during the merger with a white dwarf companion, and the single degenerate 
scenario where the white dwarf accretes material from a normal star (typically giant star) 
companion. Although simulations of double degenerate mergers suggest that some of these 



mergers produce type la supernovae (Yoon et al. 2007 Pakmor et al. 2010), it is believed 



that most of these mergers will form accretion-induced collapse implosions instead of type 



la supernovae (Nomoto & Kondo 1991; Nomoto et al. 2001 Livio 2001). Based on stellar 



accretion models, we would assume that the primary progenitor of type la supernovae is the 
single degenerate accretion scenario. 

However, population synthesis calculations currently argue against the single-degenerate 
scenario. To date, population synthesis calculations predict roughly an order of magnitude 
more double than single degenerate scenarios (see Livio 2001 for a review). The double 
degenerate merger rate is consistent with the observed supernova rate (hence, the single- 
degenerate rate is an order of magnitude too low to explain observations). 



The rate is not the only observational constraint we have on the progenitor of type la 
supernovae. The spectra from supernova outbursts provide clues as well. Hoflich (2005) has 
used the light-curve and spectral observations of supernovae to constrain the progenitor, the 
environment and the details of the explosion. Among these, one can constrain the level of 
asymmetry and the mass of the white dwarf. Both sub- and super-Chandrasekhar white 
dwarfs produce features (e.g. slow rise/decline in the super-Chandrasekhar white dwarfs) 
that differ from standard type la supernovae, arguing that most la explosions are produced in 



the explosion of a Chanrasekhar-massed white dwarf ( Mazzali et al. 2007). In this paper, we 
combine population synthesis studies of binary mergers, hydrodynamical models of mergers 
and radiation-hydrodynamics calculations to study the observational implications for the 
double degenerate Type la progenitor scenario. This combined theoretical effort allows us 
to make definitive observational predictions for the double degenerate scenario. Comparing 
these predictions with observations can place limits on the fraction of Type la supernovae 
produced by double degenerate progenitors. 

The codes used in this paper have been detailed in other papers. In Chapter 2, we briefly 
describe the codes used and outline how we couple the codes together. This coupling of codes 
allows us to make stronger predictions than have been made in the past. Most important is 
the coupling of stellar merger models, guided by population synthesis models, to our light- 
curve calculations. We describe the fate of these mergers and how they couple to light-curve 
calculations in Chapter 3. The environment created by these mergers drastically affects the 
emission from these explosions. In chapter 4 we study the effects of this environment on 
the initial shock emergence. Chapter 5 discusses the light-curves and spectra around peak 
for these explosions. We conclude with comparisons to supernova observations and their 
subsequent constrants on the double degenerate scenario. 



2. Code Description 

Our populations synthesis studies are from the simulation results of Ruiter et al. (2009) 
which used the Star Track population synthesis code. A detailed description of the algorithm 
and the input physics is given in Belczynski et al. (2002, 2008). Ruiter et al. (2009) found 
that, assuming all CO+CO mergers above a Chandrasekhar mass produced la supernovae, 
the double-degenerate scenario dominated the total, time-integrated la rate (although, de- 
pending upon the model, the rate of single degenerate mergers can be comparable after 
about 2 Gyr after a burst of star formation). Here they assumed ONe white dwarfs above 
a Chandrasekhar mass collapsed to form neutron stars and there were no CO-He mergers 
with masses above the Chandrasekhar limit. In this paper, we use the distributions of the 



total system mass from these simulations to guide our type la progenitors. 



Based on the masses in our population synthesis, we model the merger of two CO cores 
using the SNSPH code ( Fryer et al.|2006 ). This code follows the evolution of the merger from 
the onset of Roche- lobe overflow through the disruption of the companion white dwarf. These 
simulations provide us with a realistic description of the distribution of matter surrounding 
the white dwarf after the dynamical phase of the merger. Analytic estimates of the cooling 
(including convection) in this post-merger structure allows us to place limits on the density 
structure at the time of explosion. We use this range of structures to calculate the range of 
light-curves and spectra we expect from these systems. 



Our supernova is produced in a gravitational confined detonation model (Plewa et al. 



2004) from the Chicago FIASH team (Meakin et al. 2009). The explosion simulation uses 



the FLASH code (Fryxell et al. 2000) and is modeled in 3-dimensions, which we have then 



mapped into 1-dimension. A comparison of our 1-dimensional explosion velocity profile on 
top of the 3-dimensional profile is shown in figure [1} Although we are modeling the general ra- 
dial profile of this explosion with our mapping, this explosion is clearly multi-dimensional and 
any comparison of this explosion model to data will require multi-dimensional light-curves. 
However, this explosion provides a modern explosion with which to do our comparisons. The 
corresponding yields for our 1-dimensional mapping are shown in figure [2] We reduced the 
yields to 14 representative elements for our opacities: H, He, C, O, Ne, Mg, Si, S, Ar, Ca, 
Ti, Cr, Fe, and ^^Ni. The total ^^Ni yield in this explosion is O.TMq. For most of the models, 
the explosion energy is roughly 1.6 x 10^^ erg, but we included a series of weak explosions 
as well 0.4 x 10^^ erg. 

The results from the above codes are all brought together in our calculations of the 



emission using the LANL supernova light-curve code Fryer et al. (2009). These calculations 



couple a radiation-hydrodynamic simulation using the RAGE (Fryer et al. 2007 Gittings 



et al. 2008) code with post-process spectra. Both codes use the LANL OPLIB database 



(Magee et al. 1995) for opacities. Although the RAGE code is capable of running with a 



large set of photon energy groups, the current calculations use a single group Rosseland- 
mean averaged opacity using the OPLIB database. For the post-process calculations, we use 
the full frequency-dependent opacity information (14,900 groups) from this database. The 
1-dimensional structure is mapped into a multi-dimensional profile (each zone is divided into 
80 angular bins). By choosing a line-of-sight and calculating the emission along the entire 
exploding star, we model the entire ejecta including limb effects. This approach also allows 
us to calculate the full Doppler effect of the expanding material, including red- and blue- 



shifted opacities. For more details, see Fryer et al. (2009). With these calculations, we can 



compare our different merger environments to predict observational trends between these 



models. 



3. Stellar Merger and Subsequent Accretion 

In this paper, we use the results from the population synthesis models from Ruiter et 
al. (2009) to guide our initial progenitor mass distribution. Figure |3] shows the results of 
three simulations using different parameters for the common envelope evolution: model 1: 
a prescription [J where a x A = 1, model 2: a prescription where a x A = 0.5, model 3: 7 
prescription where 7=1.5. The 7 prescription of model 3 is similar to that of Nelemans & 
Tout (2005)r| the difference being that we assume this 7-formalism for common envelope 
evolution every time a common envelope event is encountered. Nelemans & Tout (2005) use 
both a 7 and a depending on the nature of the common envelope (the a-formalism is never 
used in our population synthesis model 3). 

There are a number of different formation channels producing CO+CO mergers. Chang- 
ing the common envelope efficiency and/or prescription changes both the rate and the rela- 
tive importance of these channels. Overall, our model 1 produces the most massive CO-I-CO 
mergers. Model 2, with its lower common-envelope removal efficiency, produces smaller orbits 
after the first common envelope phase. Many of these systems merge before the formation of 
the second white dwarf (see Ruiter et al. 2009 for details). These two models are standard 
models from Ruiter et al. (2009) and we will focus on their results in this paper. Although 
the two simulations predict different rates and merger times, they have several similarities. 
Both mass distributions are double peaked, one ~ I.O5M0, the other ~ 2M0. If we restrict 
ourselves to only those systems with total masses above the Chandrasekhar limit, over a 
third of the systems lie in the peak at ~ 2M0. 

Our population synthesis caluclation using the 7 prescription (model 3), however, pro- 
duces even fewer massive CO+CO mergers than our a prescription models. This is mainly 



^Thc common envelope follows a prescription: acE ( - 



GA/do„.finMp,i„„y _ GMdon.intMpHma 

''°"''°* — ''°"''""' where G is the gravitational constant, Mprimarv: -^^don int and A^don fin are the mass of 

the primary, the initial mass of the donor and final mass of the donor stars respectively, Ajnt and A^^ are 
the initial and final orbital separations, -Rdon.iob is the donor Roche lobe radius and A is a measure of the 
central concentration of the donor (Webbink 1984). 

^The common envelope is treated as a evolution of the angular momentum as mass is lost from the system: 
AJ/J = 7AM/M where AJ and AAf are the change in the angular momentum and mass respectively, J 
and M being the initial angular momentum and mass. 7 is the efficiency at which angular momentum is 
lost with mass ( |de Kool|[l990{ [Nelemans eral||2000[ [Nelemans fc Tout||2005 ). 



due to the fact that the dominant formation channel which contributes to CO+CO white 
dwarf mergers at Mtot >~ 1.8Mq in the first two models is nearly absent in model 3 (pro- 
ducing no "secondary spike" of events ~ 2M0). We will discuss this further in chapter p] 

With the peak in total system masses around 2M0, we modeled the merger of a 0.9 Mq 
WD with a 1.2 Mq WD. We follow this merger until the less-massive white dwarf is com- 
pletely disrupted. This leaves behind a compact core with an extended merger envelope 
(Fig. H). As with our explosion model, the structure is highly aspherical. Nonetheless, for 
this study, we map this structure into a 1-dimensional profile (Fig. [s]). As the core accretes 
above the Chandrasekhar limit, it will collapse (and possibly produce a type la explosion). 
Above this white dwarf is the rest of the debris from the merger. The density profile of 
this debris has a good deal of structure, but in 1-dimension, it is reasonably well-fit by: 
p = Pair Itq)'^ where r is the radius from the center and po is the density at radius tq. 

For the system modeled in our merger, the inner 1.4 Mq will accrete onto the white 
dwarf fairly quickly. By using the conditions at the end of our merger calculations, we 
estimate the cooling time (set to Erad/Lrad where Erad is the energy in this inner atmosphere, 
dominated by radiation, and Lrad is the radiative flux out of this inner material) to be 
< 10, 000 s. As this material cools, it accretes onto the white dwarf, ultimately pushing it 
above the Chandrasekhar limit. If this produces a thermonuclear explosion, it will drive a 
shock through the remaining 0.7 Mq. Due to the short cooling time of the inner material, the 
density profile of this outer 0.7 Mq will be very similar to its initial state (its cooling time is 
many orders of magnitude longer). For our light-curve calculations, we use this r~^ profile. 
With less massive mergers, the characteristics of the density profile might change slightly, 
but this r~^ profile should not change. We use this profile for the bulk of our simulations. 
If the time-scale is much longer and the atmosphere reaches an equilibrium, entropy-driven 



instabilities will produce a constant entropy atmosphere (Houck & Chevalier 1992 Colgate 



et al. 1993). For these compact object systems, the atmosphere is reasonably well described 
by an r~^ profile and this represents the shallow extreme of the possible density profiles. 
To test the dependence of our results on this choice, we model one calculation with an r^^ 
profile. 

Assuming this power law density distribution, we can then set the initial conditions of 
our models as a function of total envelope mass. For our light-curves, we use 3 different 
envelope masses: 0.1, 0.35, and 0.7 Mq. This corresponds to total binary masses of 1.5, 1.75, 
and 2.1 Mq. The companion star could be a He white dwarf or a CO white dwarf. We set 
the composition to these two separate abundances (helium vs. carbon and oxygen)r} The 



^These helium masses are not included in our population synthesis models and the higher masses are not 



entire suite of models is summarized in Table 1 



4. Shock Emergence 

Our radiation-hydrodynamics simulations are ideally suited to calculating the initial 
emergence and early evolution (through peak) of these buried type la supernovae. The 
radiation-hydrodynamics simulations allow us to include any additional shock heating as the 
exploding white dwarf interacts with its surroundings. Simulations that calculate radiative 
transfer on a homologous outflow scheme will not include this shock heating. In addition, at 
shock breakout, the material temperature is not set entirely by the radiation. This disequi- 
librium betwen radiation and matter can not be modeled by pure radiation calculations. 

We use the density/temperature profiles of these calculations to provide input for our 
detailed (> 10, 000 energy group) calculations of the spectra. With these spectra as a 
function of time, we can calculate the light-curves in any energy band (see Fryer et al. 2009 
for details). 

As the type la supernova expands, the density at the head of the shock decreases. The 
first visible emission of the supernova explosion occurs when this density becomes so low that 
radiation is no longer trapped in the flow. Hofiich & Schaefer (2009) recently studied this 
emergence of the supernova shock. For a normal type la supernovae, they argued that this 
occurs when the shock is roughly at 9000 km. At shock emergence, the type la emits a burst 
(up to 10^°ergs~^) of X-ray and gamma-ray radiation. Piro et al. (2010) argue that the 
emission peaks closer to lO'^^ergs"^. In both cases, the emission has dropped considerably 
by 1 s. 

This shock emergence is similar to shock breakout of type Ib/c and II supernovae in 
that the initial emission is driven by the thermal energy in the shock. The temperature in 
the shock tends to be high, resulting in an initial burst of high-energy photons. The primary 
difference is that the progenitor stars of type II supernovae are more extended and the 
radiation remains trapped until the shock breaks out of the star (10^ — 10^ km). Type Ib/c 
stars are more compact, but their strong winds are sufficiently dense to trap the radiation 
until the shock is further out as well. At these higher radii, the thermal temperature is much 
lower, and the breakout emission tends to be at much shorter wavelengths (UV and X-ray 
bands). The corresponding emission timescale is longer (Frey et al. 2010). 



physical. But these models give us an idea of how composition can affect the spectra and light-curves. 
^This table also includes many of the observational results discussed in chapter 5 



The type la SNe from WD mergers studied in this paper are very similar to the shock 
breakout in Type II supernovae. The merger debris acts similarly to a Ib/c wind, trapping 
the radiation until the shock is well beyond the white dwarf itself (~ 10^'^ cm for our steep, 
r~^ density profiles and > 10^^ cm for the shallower r~^ profiles). For explosions on a bare 
white dwarf, the radiation becomes untrapped when the shock breaks out of the white 



dwarf, roughly 10 cm (Hoflich & Schaefer 2009). This leads to a longer emission timescale 



and a smaller peak emission in the X-ray than that predicted for normal ("bare") Type la 
supernovae (Hoflich & Schaefer 2009; Piro et al. 2010). Figure p^ shows the total emission 



above 300eV for 5 of our models. In most cases, the peak in the X-ray flux arises within 
half a day of the explosion. The peak emission is roughly lO^^ergs"^ and this peak lasts 
10,000-30,000 s (0.11-0.35 d). As the surrounding material increases, so does the duration of 
the X-ray emission. But the structure of the density profile plays a large role in determining 
the peak flux and duration. Note that the model with 0.35 M© surrounding material with 
the r~^ dependendence on the density has a much later peak (2.8 d after explosion) and a 
peak flux of only lO^^ergs"^. 

These trends are easy to understand based on the mass of the surrounding atmosphere. 
Figure [7] shows the velocity, density and temperature structures of 4 of our models at shock 
breakout. We expect the breakout to occur when the shock hits the photosphere (or roughly 
when the mean free path is on par with the width of the shock). The width (Ar) of the 
shock grows with time, roughly staying a constant fraction of the shock size or position: 
Ar oc Tshockbrcakout- The mean free path is given by (crp)^^ where a is the opacity in cm'^ g^^. 
Setting the mean free path equal to the shock width and using our density profile (p = 
^'^shockbreakout whcrc k dcpcuds linearly on the envelope mass) , we find: 

'^shockbreakout OC '"shockbrcakout I \'^^) ■ ( 1 ) 

Solving for Tshockbrcakout , we obtain: 

''"shockbreakout OC \rx,a) . \1) 

The ratio of the shock breakout radii (rshockbrcakout.Mi/'^shockbreakout.Ma) for two different enve- 
lope masses (Ml, M2) is simply: 

^shockbreakout,Mi/''^shockbreakout,M2 = (^Ml/^M2) = (M1/M2) . (3) 

We expect, then, that the shock in the 0.7 Mq model to experience shock breakout roughly 
1.9 times further out than the 0.1 M© model. The duration of the breakout should be roughly 
proportional to r/vr where Vr is the shock velocity. Hence, we also expect longer durations 
for our more massive atmospheres (From figure [TI we see that both the radius is bigger and 
the radial velocity is smaller). 



The primary outlier in this picture is the shallower density profile model. But the same 
physics explanations apply. The shallower density means that shock breakout occurs at 
much higher radii. Not only is the breakout delayed (and the duration of the breakout is 
longer), but at these large radii, the shock has cooled, leading to a much weaker signal. To 
accurately calculate the breakout emission, we will need to have accurate estimates of the 
density profile of the surrounding material. 

Can we use shock breakout emission to constrain the progenitor? If the breakout emis- 
sion is extended for more than 1000 s, we are assured that non-neglibile material lies on top 
of the exploding white dwarf. Technically, the duration scales with the amount of mass on 
top. But as we have found from altering our density profile, the profile is as important, if 
not more important, than the total amount of mass. Until we have accurate models of the 
merger process and can model the explosions through these profiles in multi-dimensions, it 
will be difficult to constrain the exact mass of the system from shock breakout alone. 



5. Spectra and Light-Curves Near Peak 

We are not limited solely to the shock breakout light-curve. The surrounding material 
will also affect the light-curve and spectra near peak. Especially for these enshrouded type 
la, shock heating remains an important feature of the light-curve through peak and our 
radiation-hydrodynamics calculations are ideally suited for calculating spectra and light- 
curves of these supernovae. 

Our basic suite of light-curve models consists of 3 different envelope masses: 0.1, 0.35, 
and 0.7 Mq, corresponding to total binary masses of 1.5, 1.75, and 2.1 Mq. For each envelope 
mass, we model two simulations, one assuming the envelope is composed of carbon and 
oxygen, the other assuming it is all helium (see Table fllfor details). In this way, we study not 
only the double degenerate progenitor consisting of two CO cores, but also the dependence 
on the composition of the debris. Figure |8] shows the spectrum over a broad wavelength and 
flux range for each of our 6 models. The corresponding light-curves for a range of wavelength 
filters are shown in Figure |9| 

Let us review the trends in the light-curves first. One pervasive feature of a surround- 
ing environment is that it can delay initial photon emergence. The surroundings trap the 
photons longer, delaying the UV burst. In addition, shock heating as the explosion hits 
this surrounding environment produces high temperatures and longer (or second) UV out- 
bursts. This means that the UV outburst may still be strong at peak V-band emission. A 
large fiux shorter than 1500 A at peak V-band is a strong indicator of a large surrounding 
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environment and, hence, a double-degenerate progenitor. If we further increase the mass 
of the surrounding medium, the radiation is trapped longer. This will produce a broader 
light-curve. But shock heating plays a larger role in driving the supernova, and this extra 
energy source makes it more difficult to analyze the light-curve. 

Many of these trends can be seen in figure [9j With a 0.1 Mq surrounding environment, 
the light-curves of the helium and carbon/oxygen mergers are very similar. The V-band 
peaks between 15-30 d with an absolute magnitude dimmer than -17. The UV emission is 
initially bright at first emergence from the expanding ejecta and then peak again due to 
shock heating (20-40 d). The peak UV emission predicts absolute magnitudes near -17 in all 
the Swift bands. 

As the environment mass increases, so does this peak emission. With a 0.35 Mq sur- 
rounding environment, the UV emission peaks at absolute magnitudes of -20. Shock heating 
is dominating this UV emission. For a helium envelope, where the opacity is lower because 
the helium in the envelope is quickly ionized, the UV emission is broadened and peaks at 
-20. The V-band also peaks higher (-19), but decays in 10 d by over a magnitude. For C/0 
envelopes with their higher opacities, the UV peaks high (-20) at later times than its helium 
counterpart. In addition, the V-band is strongly broadened. Even after 60 d, it is within 0.2 
magnitudes of its peak. 

At still higher envelope masses (0.7 Mq), photons are trapped longer. When photons 
do escape, the temperature of the ejecta is lower, leading to lower UV emission than in 
our 0.35 Mq case. Even for helium surroundings, the V-band emission begins to produce a 
broader peak. But this broadening is exacerbated in the high-opacity C/0 surroundings. 
For massive C/0 envelopes, the V-band is still rising after 60 d. 

The behavior of these light-curves is very different from those of normal Type la super- 
novae. These differences, in part, can be understood by the differences in the structure of the 



star throughout peak. Figure [TO] shows the velocity, density and temperature profiles of our 
standard 0.35 Mq and 0.7 Mq envelopes at a range of times during the peak breakout emis- 
sion. Note that the shock remains at the boundary between trapped fiow and free-streaming 
throughout peak. The front of the shock is pushing through low-density media and emitting 
nearly at the free-streaming limit. But the width of the shock as defined by the temperature 
profile is roughly a mean-free path thick. This is very different than an inner, optically thick 
blackbody source assumed in bare type la supernovae. In addition, although the velocity 
profile approaches a homologous outfiow, it is not quite a homologous outfiow. Small varia- 
tions in the shock temperature can cause large variations in the emission, especially the UV 
emission. 
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This complexity, especially the fact that shock heating is playing a role in the light-curve, 
makes it difficult to make detailed predictions from enshrouded type la supernovae. There are 
a few robust claims. First, enshrouded la supernova (with more than 0.1 Mq of surrounding 
material) should have broader V-band light-curves. The light-curves are, in part, powered 
by shock heating, and these shocks lead to higher UV emission. Thus far, the current set of 



Swift UV type la supernova light-curves do not show this high UV emission (Brown et al. 



2009), suggesting that none of these observed outbursts are enshrouded type la supernovae. 
Strong UV emission and/or extremely broadened V-band peaks are tell-tale signatures 



of massive C/0 surroundings. Figure 11 shows the full spectra for our models. Near peak, 
we see the strong calcium and silicon lines we expect in a type la supernova. With the 
uncertainties in the explosion and the outer density structure, it is possible to produce 
spectra at peak that ffi many observed type la supernovae. 

But the spectra at early times are much more sensitive to the surrounding medium. For 
example, a broad absorption feature at ~ 3700A due to Call H&K and Si II is present at 



early times in observed type la supernovae Bufano et al. (2009), but it is not present before 



peak in some of our models. The surrounding C/0 material dominates the absorption 
features near peak, lessening the absorption features by synthesized elements until peak. 
Early observations may be the key to finding distinguishing spectral features. 

We have run a series of other models, varying the explosion energy and the density 
profile. These models all show the same trends, extended light-curves and enhanced UV 



emission. For example, figure 12 shows the light-curves and spectra of two of these models: 
a weaker explosion and a shallower density profile. The peak in the weak explosion occurs 
later than our standard energy explosion, but the duration of the peak (e.g. the time the 
absolute visible magnitude is above -18) is roughly comparable to the standard explosion. 
A shallower density profile produces a dimmer, but much longer duration, outburst. 

We end with one cautionary note on exact line comparisons. A major difference in 
our models is that our post-process spectral calculations are built on top of a true "two- 
temperature" radiation hydrodynamics calculation. By "two-temperature", we refer to the 
fact that the material temperature is affected by shocks as well as the radiation. In typical 
type la supernova calculations, shock heating is neglected (it is argued that does not play a 
role in the observed emission for photons with wavelengths longer than the X-ray or U V) . For 
these enshrouded progenitors, shock heating is more important and our models include this 
effect. But our models currently lack the detailed resolution at the photosphere to catch line 
features sensitive to the narrow region just beyond. And we have not yet included the fact 
that the excitation levels of each atom are not described by a single equilibrium temperature. 
Although the broad features discussed here will not change, many of the details of the spectra 
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will change as more accurate models are made. 



6. Putting it all together 

In this paper, we have combined population synthesis studies of the double degenerate 
scenario for type la supernovae with merger calculations of these CO/CO binaries to produce 
the environments surrounding this type of type la supernova. We have then modeled the 
propagation of a modern explosion simulation through this environment to produce light- 
curves of these enshrouded type la supernovae. The surrounding environment both delays 
and extends the X-ray flux in shock breakout, producing an initial X-ray signature much 
closer to the shock breakout of a type Ib/c supernova than that seen in the shock emergence 
of a bare type la. The subsequent V-band peak is extended and is much broader than typical 
las. Finally, although the peak and post-peak spectra will display the same lines as normal 
supernova, the early spectra will be dominated by carbon and oxygen lines only. 

We have made a number of approximations in our calculations. We assumed that the 
density profile maintained its steep, r~'^ profile. If the surrounding atmosphere redistributes, 
we would expect a slightly flatter profile. Our test model with a r~^ profile produced an initial 
X-ray flash that occurs later and the visible light-curve is broader and dimmer. Flattening 
the density proflle produces an explosion very different than typical la supernovae. 

Another approximation we made was to spherize the inherently multi-dimensional merger 
and explosion models. The merger debris is asymmetric, slightly denser in the orbital plane. 
The light curve will be affected if the explosion is also asymmetric, sending a weaker explosion 
along the orbital plane. To test this affect, we modeled a weaker explosion. Although the 
light-curve peaked later, it was similarly broad to our strong explosion, again not matching 
the typical type la supernova. 

All of these results argue that these enshrouded systems can not explain normal type 
la supernovae. And, such systems must be less common than normal type la supernovae 
to avoid dominating our current la supernova sample. Pakmor et al. (2010) argued that 
mergers, with a weak explosion, might explain sub-luminous mergers. They found that the 
only issue with such an explanation is that the merger model produced light-curves that were 
slightly broader than the observations. If they had run radiation-hydrodynamics calculations 
(including the energy from shock heating), it is likely that the light curves would be even 
broader, making an even worse flt to sub-luminous models. 

Astronomers have argued for other forms of enshrouded type la systems, especially 
for the case of Supernova las with hydrogen envelopes. For example, supernova 2002ic is 
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fit well by models invoking a wind profile with mass loss rates in excess of 0.01 MqY ^ 
or circumstellar material extending out to 10^^ — 10^^ cm with a possible gap between the 



exploding white dwarf and this circumstellar material ( Chugai et al. 2004 Deng et al. 2004 



Kotak et al. 2004 Wood-Vasey et al. 2004). Supernova 2005gj may be similar (Aldering et 



al. 2006). Our surrounding debris is much more compact, concentrated below 10 cm. 



Some systems exist with strong carbon lines, suggesting incomplete C/0 burning (Maz 



zali et al. 2001; Tanaka et al. 2008). Typically these systems have considerable silicon from 



this partial burning. If the envelopes of double degenerate mergers are particularly compact, 
we might be able to match these systems, but we haven't studied such systems here. 

But if double degenerate systems are to dominate the la sample, why don't we observe 
them? 

One explanation for the lack of observations of these enshrouded type la supernovae is 
simply that CO/CO mergers do not produce type la supernovae. Nomoto and collabora- 



tors (Nomoto & Kondo 1991; Nomoto et al. 2001) have argued that the high accretion rate 



in these systems will preferentially lead to the collapse of the merged white dwarf system 
to a neutron star. The resultant explosion will still produce a supernova-like explosion, but 



could be much dimmer than what we studied here (Fryer et al. 2009) 



However, other possible explanations exist. Our current stellar models predict too-large 
CO white dwarf masses. There is some evidence that the ^^C/^^C fusion rate has been under- 
estimated. By increasing this rate, Herwig et al. (in preparation) found that the maximum 
CO white dwarf mass is much lower (as low as 0.8Mq). In such a case, the maximum total 
merger mass from CO binaries would be below 1.6Mq. A number of ONe binaries would 
exist in this scenario, but these mergers are expected to collapse to neutron stars, yielding 
the dimmer supernova discussed in the previous paragraph. 

Another potential explanation is that our standard population synthesis models predict 
a peak at around 2M0 in the total mass distribution of type-la producing CO/CO mergers, 
arguing for massive surrounding envelopes (~ 0.6Mq). But by modifying the prescription 
for binary mass transfer, i.e., by adopting the 7-formalism (model 3) rather than the a- 
formalism, the peak at 2 Mq disappears, arguing that most mergers with masses above 
the Chandrasekhar limit occur with masses within 0.1 Mq of this limit [3j Without a better 
understanding of common envelope evolution, population synthesis models may simply be 
unable to predict the mass distribution of merger WD binaries accurately. Perhaps there are 
simply few enshrouded type las and most systems only have 0.1 Mq of surrounding material 
or less. If such systems can explode to form type la supernovae, their observations would 
not be so distinct from the current type la population. 
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What we can say for certain is that the spectral and hght-curve features of double degen- 
erate mergers with total system masses in excess of roughly I.S-I.GMq do not match normal 
type la supernovae. These super-Chandrasekhar mergers can not dominate la progenitors. 

This work was carried out in part under the auspices of the National Nuclear Security 
Administration of the U.S. Department of Energy at Los Alamos National Laboratory and 
supported by Contract No. DE-AC52-06NA25396. 
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Table 1. Transient Models 



Name 


MEnvelope 


Envelope 


Vpcak 


Vfwhm 


UVpeak 


uv 




M© 


Composition 


Abs. Mag. 


d 


Abs. Mag. 
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Double Peak 
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0.35CO 
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0.35he 
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He 
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-20 


Broad (40 d) 


0.7co 


0.7 


CO 


-17 
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Delayed (40 d rise) 


0.7he 


0.7 


He 


-17 
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Broad (>40d) 


0.35CO Weak 


0.35 


CO 


-19 


40 


-20 


Delayed (60 d rise) 


O.lco Weak 


0.1 


CO 


-17 


40 


-17 


Delayed (100 d rise) 


0.35co-r3 


0.35 


CO 


-19 


> 100 


-20 


Broad 




109 

Radius (cm) 



lO'o 



Fig. 1. — Velocity profile of our type la supernova explosion, based on simulations by Meakin 
et al. (2009). The black points depict the velocity versus radius from this simulation. The 
line shows our 1-dimensional fit to these data. Our 1-dimensional fit can not match the full 
structure of this explosion, but it does match the basic structure and conserves the total 
energy. 
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Fig. 2. — 1- dimensional abundance profile of our type la supernova explosion (see Fig. IT] 
for details). We have combined the yields to focus on the major trends in this distribution 
profile. For our opacities we use 14 representative elements. 
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Fig. 3. — Distribution of total masses for our population of merger binaries based on the 
population synthesis models of Ruiter et al. (2009). We will assume only those binaries 
whose total mass exceeds 1.4 M© will produce explosions. Note that the peak in masses 
above 1.4 Mq occurs at roughly 2M0. 
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Fig. 4. — A slice in the x-y (orbital) and x-z (out of orbital) plane of our 3-diniensional 
simulations of white dwarf mergers (a simulation using the SNSPH code). The coloring 
denotes density (blue is high, red is low) and the vectors denote velocity magnitude and 
direction. Although far from symmetric, we will angle average this density profile for our 1- 
dimensional explosions. Clearly, multi-dimensional models are required to produce accurate 
density profiles and, presumably, spectra and light-curves for the supernova in these systems. 
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Fig. 5. — The density profile of tlie outer envelope of our star. Although there is considerable 
structure in this profile, it is fit reasonably well with an r~^ density structure. We will assume 
this power-law for most of our simulations. 
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Fig. 6. — The X-ray luminosities as a function of time from shock break-out for 5 of our CO 
surrounding atmospheres: O.lCO-Weak (sohd), 0.35CO-Weak (dotted), 0.35CO (dashed), 
0.35CO-r3 (long-dashed), 0.7CO (dot-dashed). Most of our models peak around lO^'^ergs"^ 
and the emission lasts for 0.1-0.3 days. The exception is the atmosphere where we used the 
r~^ density profile. In this model, the radiation remains trapped longer and shock break-out 
occurs at higher radii when the shock is cooler. This produces a lower X-ray flux with a 
longer duration. All of these models are very different from the shock emergence seen in 



normal la supernovae, which have no surrounding atmospheres (Hofiich & Schaefer 2009 



Piroet al. 2010). 
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Fig. 7. — The velocity (top), density (middle), and temperature (bottom) profiles at shock 
break-out of 4 of our shock models: 0. ICO- Weak (solid), 0.35CO-Weak (dotted), 0.35CO-r3 
(long-dashed), 0.7CO (dot-dashed). Shock breakout occurs later for more massive envelopes. 
The shallower density gradient for model 0.35CO-r3 traps the radiation longer. Shock break- 
out for this model occurs after the shock has traveled 15-30 times further out. This explains 
the later X-ray emergence and shallower decay in the X-ray luminosity. Its temperature is 
correspondingly lower at shock-breakout, leading to a lower X-ray luminosity. 
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Fig. 8. — Spectra (luminosity versus wavelength) from our 6 primary explosion models at 
a series of time positions. The left column corresponds to helium atmospheres, the right 
column corresponds to CO atmospheres. The envelope masses range from 0.1 to 0.35 to 
0.7 M0 for the top, middle, and bottom panels respectively. At early times, there is a peak 
at low wavelengths (high energies) denoting the initial emergence of the radiation. With time, 
this peak moves to lower energies (longer wavelengths) . Ultimately, the material becomes so 
diffuse that it is ionized by the radiation and many of the red line features disappear. 
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Fig. 9. — Light-curves (magnitude versus time) of our 6 primary models for a range of filter 
bands (See Fig. [s] for details). Here we are using the filters used in the Swift satellite. Shock 
heating in the envelope produces extended UV emission and may also broaden the emission 
in the V-band. As the envelope mass increases, the emission is first increased (as shock 
heating becomes more important), and then decreases as photon trapping limits and delays 
the emission. 
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Fig. 10. — The velocity (top), density (middle), and temperature (bottom) profiles near shock 
break-out at a variety of times for our more massive CO envelope models. Throughout the 
peak emission for both models, the shock remains at the boundary between trapped radiation 
and free-streaming radiation. 
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Fig. 11. — Spectra for our 6 primary models (See Fig. p^for details) focusing on the optical 
and infra-red spectra. 
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Fig. 12. — Light-curves (magnitude versus time) and spectra at specific times for 2 of our 
additional models: 0.35CO-Weak and 0.35CO-r3. Note that these models also have the long 
delays and enhanced UV emission. 



